Network Neuroscience

Micro-electrode array (MEA) recordings are widely used to characterize functional connectivity in neural cultures and have gained traction for the analysis of human brain slices. However, the impact of graph construction methodology on the resulting network topology has not been systematically quantified. Here, we benchmark three methods - shared spiking activity, Pearson cross-correlation, and the spike time tiling coefficient (STTC) - across 37 recordings from human cortical slice cultures classified into low, moderate, and high activity groups. We show that method choice alone produces large topological differences (Cohens d = 0.86-1.14 for clustering coefficient, d > 1.0 for node count), while higher-order features such as modularity remain stable. Each method exhibits a distinct sensitivity profile: shared spiking detects activity-dependent changes primarily through network size, correlation uniquely captures clustering differences, and STTC combines strong biological sensitivity with negligible parameter dependence across lag windows (all d < 0.1). Within shared spiking, z-score normalization dominates all other parameter choices (d > 1.0 versus bin size effects of d < 0.23), functioning as an implicit analytical null model that fundamentally reshapes the edge set rather than merely rescaling weights. Inter-method edge overlap is low (Jaccard index 0.08-0.45) and activity dependent, demonstrating that these methods identify substantially different connections from identical data. Our results reveal that methodological choices including construction method, threshold, and normalization introduce hidden degrees of freedom with effect sizes comparable to the biological signals being measured. We provide practical recommendations for parameter selection, reporting, and cross-method validation in MEA-based network neuroscience. Author SummaryWhen we record electrical activity from brain tissue using grids of electrodes, we can ask how different sites influence one another and map the tissue as a network of connections. Thanks to novel culturing methods, this approach is increasingly used to study human brain slices. However, deciding what is "connected" is not well defined. Researchers use several different methods, and it has never been clear how much this choice shapes the network they end up describing. Here we compared three widely used methods on 37 recordings from human cortical slices spanning a range of activity levels. We found that the method alone can change the apparent structure of the network as much as real biological differences do. The methods frequently disagreed about which connections exist and some technical choices, including normalization techniques, had surprisingly large effects. Because these hidden choices can rival the biological signal, we provide this benchmarking work with practical recommendations for selecting, reporting, and cross-checking methods, so that network studies of brain tissue become more transparent, comparable, and reproducible.

Brain activity can be understood as a sequence of neuronal avalanches, i.e., transient episodes of coordinated activation that emerge across scales, from individual neurons and local networks to whole-brain dynamics. Avalanches are typically characterized by features such as size, duration, number of active components, and the silent time separating consecutive events. Although these features have been extensively characterized through their marginal distributions, their temporal organization and dependence on the underlying brain architecture remain poorly understood, leaving us without a framework for embedding fast neuronal avalanches within slower brain dynamics. Here, we analyzed eyes-closed resting-state magnetoencephalography recordings and the corresponding structural connectomes from 30 healthy participants to investigate the dynamics of avalanche sizes. We found that large avalanches preferentially followed short silent times, whereas small avalanches were more likely to occur after long silent periods. Based on the empirical joint distributions of avalanche size and silent time, we could define four types of events occurring above chance levels (avalanche large or small, preceding pause long or short). Mixed categories--combining a small value of one feature with a large value of the other--occurred more frequently than expected, while same-category events happened less often than chance. Furthermore, consecutive events tended to remain in the same category, a phenomenon referred to as persistence. We next investigated whether a brain regions connectivity profile shapes its propensity to participate in avalanches of different sizes. More strongly connected regions participated most often in small avalanches, whereas weakly connected regions were preferentially recruited during large avalanches. This pattern may reflect the greater sensitivity of highly connected hubs to fluctuations propagating through the network, resulting in frequent but spatially contained events. By contrast, the recruitment of more peripheral regions may require broader and stronger collective activity, occurring only during rarer, large-scale avalanches. In contrast, regional participation showed no clear association with the silent time preceding an avalanche. Together, these findings show that neuronal avalanches are neither temporally independent nor anatomically unconstrained: their sequence retains a memory of preceding events, while structural topology shapes which regions are recruited as avalanches grow. By connecting fast avalanche dynamics with slower temporal organization and the structural connectome, our results provide a multiscale framework for understanding how transient events are embedded within ongoing brain activity.

Objective: Brain network controllability provides a framework for understanding how structural organization shapes brain dynamics, yet current models mainly rely on white-matter connectivity and may overlook the contribution of gray-matter architecture. Approach: We constructed a fusion network combining diffusion tensor imaging-derived white-matter connectivity with gray-matter morphological similarity and investigated its controllability, biological associations, heritability, phenotype prediction, and control energy. Main results: Controllability derived from the fusion network preserved key topological properties of the white-matter network and was associated with neurotransmitter systems and cerebral metabolism. Compared with the white-matter connectivity-based network, fusion-based controllability showed a systematic shift toward higher heritability, improved prediction of several individual characteristics and cognitive functions, and lower modeled control energy for activating resting-state networks. Significance: These findings suggest that incorporating gray-matter morphological information into a DTI-supported network provides a complementary structural representation for studying brain network controllability and state transitions. The lower control energy represents a model-derived transition cost and should not be interpreted as a direct measure of physiological energy expenditure.

Understanding dynamic mechanisms underlying cognition remains a major challenge in human neuroscience. Here, we develop, validate, and apply Multivariate Dynamical Systems Identification with Amortized Variational Inference (MDSI-AVI), a novel computational framework designed to address critical challenges in capturing asymmetric, context-dependent, whole-brain directed interactions while accounting for regional hemodynamic response variability in fMRI data. MDSI-AVI leverages simulation-based inference through forward and reverse variational inference to address the limitations of conventional variational methods in high-dimensional settings. By averaging over uncertainty in hemodynamic response parameters using forward simulation, MDSI-AVI provides well-calibrated posteriors of directed connectivity that scale efficiently to networks with hundreds of nodes. Applied to Human Connectome Project data (N=728), MDSI-AVI reveals new insights into working memory mechanisms, identifying the dorsal anterior insula as a critical hub influencing activity at the whole-brain level. We demonstrate task-dependent modulation of causal influences, where the salience network drives frontoparietal network activity, which differentially influences the default mode and sensorimotor networks depending on working memory load. These whole-brain causal interactions distinguish task conditions with high accuracy and predict working memory performance. Our framework demonstrates reproducible results across whole-brain parcellations, establishing MDSI-AVI as a robust tool for advancing our understanding of circuit dynamics in cognition and disease.

Intrinsic neural timescale (INT) quantifies the persistence of spontaneous neural dynamics and offers a principled metric for characterizing brain-wide temporal organization. Although a hierarchy of INTs has been established during rest, how task engagement reconfigures this organization and how it is constrained by the structural connectome (SC) remain poorly understood. Here, we systematically mapped whole-brain INT using high-resolution fMRI data from the Human Connectome Project during rest and seven tasks spanning working memory, gambling, motor, language, social, relational, and emotion domains. Task engagement induced robust, regionally heterogeneous changes in INT while largely preserving the brain-wide temporal hierarchy across cognitive states. SC-INT coupling remained strong but consistently decreased during tasks, indicating that anatomical architecture continues to constrain INT, although its influence is attenuated under task demands. To investigate these findings mechanistically, we employed a multiscale, whole-brain neuronal-network model, which revealed that INT increase and peak within a broad critical-like regime. Strong SC-INT coupling, as observed empirically, emerged in the subcritical regime, weakened progressively with increasing network excitability, and reversed in the supercritical regime. These results demonstrate that task engagement reconfigures INTs while maintaining their hierarchical organization, suggesting that both resting and task states operate largely within a common subcritical dynamical regime.

Two decades of functional connectivity research have established schizophrenia as a disorder of distributed dysconnectivity, with a robust thalamocortical signature: reduced prefrontal and increased sensory coupling. A fundamental issue is that functional connectivity is undirected, operates at a single slow timescale, and cannot reveal causal direction. Moreover, the mismatch between BOLD sampling speed and neural dynamics can hide edges, fabricate spurious ones, and reverse the apparent orientation of causal relationships. To overcome these limitations, we introduce a general framework for estimating directed causal graphs from fMRI that explicitly accounts for temporal undersampling. We apply RnR, a causal discovery method built on the rate-agnostic RASL framework, to resting-state fMRI from the multi-site FBIRN cohort. Rather than returning a single directed graph, RnR recovers an equivalence class of directed graphs consistent with the observed data, each annotated with the sampling rate that would produce it and classifies each estimated orientation by its stability across inferred rates. This provides a principled basis for distinguishing directed interpretations that are safe to trust from those that are timescale-contingent. Benchmarking against five single-timescale estimators on schizophrenia data, RnR recovered substantially more group-differentiating directed edges, reproducing the fields most replicated finding in directed form: a sensory-to-visual hyperconnectivity hub oriented from the post-central gyrus component to primary visual cortex. Through simulation, we show that coarse spatial resolution has shielded single-timescale methods from the full effects of undersampling, whereas finer parcellations will require explicit undersampling modeling. This work reframes fMRI effective connectivity estimation by treating undersampling as a fundamental property of the measurement, enabling directed interpretations that are grounded in the data-generating process.

Network measures of the ageing connectome are dominated by magnitude: connection strength and density decline, and the topological summaries built on them decline with them. Whether the geometry of the network follows the same course is not known, because the quantities in common use do not separate how strong a connection is from how it sits among the connections around it. We computed the Ollivier-Ricci curvature of every edge in structural connectomes from 307 participants spanning the adult lifespan, a quantity defined by optimal transport between the neighbourhoods of connected regions, and asked how it changes with age. The total geometric separation between within-network and between-network connections did not change across seven decades. Underneath that constancy, individual network pairs moved substantially and in opposite directions, gaining curvature around the salience and ventral attention system and losing it between the control and default mode networks. Curvature and connection strength reached half of their age-related variation almost four decades apart, and a small set of prefrontal nodes moved against the global trend. Ageing appears to conserve the local redundancy of the structural connectome in total while relocating it, on a timescale distinct from that of connection strength.

The gastric network, comprised of brain regions whose activity synchronizes with the stomach's slow-wave rhythm, offers a unique window into the brain-body interaction involved in interoceptive processing. While previous work has established the existence of this network, its intrinsic organization and temporal unfolding remain poorly understood. Here, we reanalyzed resting-state fMRI-electrogastrogram data from 43 healthy adults of both sexes to characterize the time-averaged architecture and time-varying reconfiguration of the gastric network. We identified regions exhibiting phase-locked synchronization with the stomach slow electrical rhythm (0.05 Hz) and characterized cortical parcels comprising this network. Time-averaged graph-theoretical analysis revealed a fixed unimodal organization of functional communities, with primary visual, default mode network (DMN) and dorsal attention regions emerging as the principal time-averaged hubs. Next, we applied edge-centric functional connectivity (eFC) to capture the network state during transient high-amplitude "bursts". Time-varying community detection revealed communities whose compositions formed integrative combinations of DMN, visual, attentional and control elements. Edge-derived hubs shifted away from primary visual dominancy in the time-averaged analysis, and were instead directed by DMN regions, suggesting that moments of heightened connectivity in the network are coordinated by multisensory integration rather than passive sensory processing. These findings demonstrate that the gastric network is not merely a time-averaged, sensory-bound system, but rather a flexible and dynamically reconfiguring interoceptive network whose organization is selectively coordinated by transient cofluctuation events. This work provides a comprehensive network analysis of gastric-brain coupling and reveals a temporally structured mode of interoceptive integration that may support adaptive physiological and cognitive regulation.

How the brains physical geometry gives rise to its flexible functional repertoire remains a central question in neuroscience. Here, we trained three classes of recurrent neural networks (RNNs) on a working-memory task, forming a graded hierarchy of spatial constraints: Vanilla RNNs (no spatial constraints), Masked RNNs (projection constraints limiting where information enters and leaves the network), and biophysical RNNs (bioRNNs; projection constraints and spatial embedding of the networks connectivity using the brains inter-regional Euclidean geometry). We assessed how well each RNN class predicted empirical fMRI activity without exposing them to it during training. Our results showed that bioRNNs were the only networks to successfully predict empirical brain activity and to organize their dynamics into a spatial pattern that recapitulated the brains principal hierarchy (the sensorimotor-association axis). Additionally, bioRNNs ability to predict empirical brain activity emerged along a trajectory in which geometry was laid down first, then partly traded back as the task was mastered. Importantly, brain-like topological features emerged in bioRNNs as they increased their task proficiency while maintaining their ability to predict brain activity. Taken together, our results indicate that physical geometry and cognitive inputs play distinct, complementary roles: while geometry constrains the space of possible brain dynamics, cognitive inputs determine which dynamics are expressed. They also situate topology as the scaffold through which the physically embedded brain reconciles wiring costs and computational demands.

Despite the growing use of functional near-infrared spectroscopy (fNIRS) hyperscanning to record brain activity simultaneously from interacting individuals in naturalistic settings, most analyses quantify functional connectivity separately for each channel pair. The resulting collection of pairwise estimates is difficult to integrate into a network-level characterization of intra- and inter-brain organization. Here, we present an open, configuration-driven Python toolkit that transforms preprocessed fNIRS hyperscanning time series into functional connectivity graphs. The toolkit constructs a bipartite inter-brain network for each dyad and separate intra-brain networks for each participant, computes node- and graph-level measures, and exports adjacency matrices, edge lists, analysis-ready summary tables, reproducibility metadata, and standardized visualizations. Dataset-specific parameters, including directory structure, participant naming, channel selection, epoch extraction, and edge-retention criteria, are defined in a human-readable YAML configuration file, enabling the same workflow to accommodate differently organized datasets without changes to the source code. We illustrate the pipeline using a representative recording from a mother-infant fNIRS hyperscanning dataset and present the resulting network outputs. The toolkit provides a reproducible framework for moving from pairwise functional connectivity estimates to network-level analyses of dyadic and individual brain organization.

Quantifying longitudinal changes in an individual's brain is central to the development of personalised neural biomarkers in neurology and psychiatry. However, existing approaches for characterising individual neurophysiological signatures focus on discrimination between people rather than the quantification of within-subject change. To address this, we introduce the Brain Stability Index (BSI), a whole-brain metric that quantifies the similarity between two longitudinal neurophysiological scans in a low-dimensional latent space, with reference to a normative magnetoencephalography (MEG) database. Using 276 open resting-state MEG datasets and matched synthetic data, we first characterise how finite test-retest reliability sets a noise floor on the BSI. We then demonstrate that the BSI is sensitive to graded changes in whole-brain neural change that extend beyond measurement variability. Finally, we show that Factor Analysis, by separating shared structure from feature-specific noise, makes the BSI more robust to measurement artefacts. Together, these findings establish the BSI as a robust, bounded measure of neural stability that is well suited to longitudinal monitoring in neurology and psychiatry.

Lesion Network Mapping (LNM) is a framework used for identifying symptom-related brain circuits by projecting lesion locations onto a normative connectome. Recent methodological investigations have raised concerns about the biological interpretation and specificity of the circuits derived using this method, with published LNM maps often showing high similarity across clinically unrelated conditions. Specificity testing has subsequently been put forward as the decisive step to ensure specificity to the symptom in question, accompanied by the argument that this step was not evaluated in the original methodological investigation. Yet, sensitivity testing, specificity testing, case-control LNM, permutation of group labels, and symptom-based LNM involve related operations on connectivity matrix C. We expand on specificity testing in LNM, clarify its relationship to other LNM steps and variants, and examine the persistent repetition among LNM specificity networks across studies. These considerations advance our understanding of the disease-specificity limitation of LNM and encourage the development of new methodological approaches for identifying brain circuits underlying psychiatric and neurological disorders.

Previous work suggests that different groups of neurons in orbitofrontal cortex (OFC) constitute the building blocks of a circuit in which economic decisions are formed. Here we used network inference analysis (Ising model) to shed light on the internal organization of this circuit. We examined populations of neurons recorded simultaneously, and inferred the functional couplings. We then computed a reduced, effective network (EN) where each node corresponded to an encoded variable. The EN had a recognizable structure, with enhanced couplings between input and output neurons supporting the same decision, and enhanced couplings between neurons encoding value variables with the same sign. This structure was highly reproducible across individuals and hemispheres. Importantly, it depended on the internal state of the animal and the behavioral conditions. The EN reproducibility decreased with the distance between cells but it increased with the number of cell pairs, suggesting that OFC operates as a single distributed assembly.

The development of neural connections in the brain results from a complex interplay between biological processes and mechanical forces. A key question in neuroscience is how physical forces and the mechanical properties of brain tissue influence the formation of structural connections. Here, we demonstrate that mechanical forces play an essential role in shaping the emergence of short-range connections, particularly U-shaped fibers that link neighboring regions of the cortex. Using a computational model that incorporates our "stress-dependent axon reorientation" hypothesis, we simulate how growing axons respond to the mechanical stress field generated by cortical folding. Our results suggest that axonal growth and reorientation may be strongly influenced by local mechanical cues, helping establish the organization of these short-range pathways. Supported by in vivo diffusion tensor imaging and histological observations, our findings provide a physical explanation for why these fibers predominantly adopt U-shaped trajectories, and why connections between gyri (ridges) are more prevalent than those between sulci (valleys) or spanning gyri and sulci. These results suggest that understanding the mechanics of brain folding is critical for fully explaining the formation of brain connectivity and its variations in health and disorder. Teaser: Mechanical forces during cortical folding guide the formation of short association fibers in the brain.

Individualized resting-state functional magnetic resonance imaging (rs-fMRI) is increasingly used to guide neuromodulation target selection. However, clinical scans are often short and noisy, and standard pipelines for functional network identification do not provide information about confidence of network assignment. With limited data, unstable network assignments can misdirect stimulation toward off-target regions, making it critical to know which assignments can be trusted. We developed Precision Confidence Mapping (PCM), a bootstrap-based framework that makes this uncertainty explicit and actionable. PCM repeatedly resamples the time series and reruns network detection to estimate how consistently each vertex is assigned to a given network. The resulting confidence maps can be thresholded to exclude less stable regions. We evaluated PCM across scan durations from 5 to 70 minutes using positive predictive value (PPV) as the primary measure of network-assignment precision. PPV quantified the proportion of vertices assigned to a network that received the same label in an independent within-subject 70 minute reference map. Confidence thresholding markedly improved PPV across functional networks, with the largest gains for short scan durations. Compared with standard network assignment, PCM significantly increased agreement with this independent reference. Within-subject agreement remained greater than between-subject agreement, indicating that thresholding preserved individual-specific network topography. These precision gains came with modest reductions in reference-network coverage, particularly at shorter scan durations. This tradeoff may be acceptable for neuromodulation applications that prioritize minimizing off-network assignments. By adding a reliability layer to individualized mapping, PCM supports more cautious and precise neuromodulation targeting under real-world clinical scan constraints.

Determining the amount of information a neuron receives per unit of time is key to understanding brain connectivity and how neural networks encode and transmit information. In particular, estimating this information flow by distinguishing between excitatory and inhibitory synaptic contributions is critical to understanding neural network function, as maintaining the excitation-inhibition (E/I) balance regulates neuronal excitability and circuit stability, whereas its disruption can lead to a plethora of brain disorders, including neurodegenerative and psychiatric conditions. However, because synaptic conductances cannot be measured directly, inverse methods are required to infer them from the membrane potential --a readily measurable quantity. Although partial solutions have been proposed, accurately estimating these conductances remains a significant challenge due to the complexity and diversity of the inputs. This is particularly true in the spiking regime, where neurons actively fire. In this work, we introduce a novel computational strategy that combines two critical metrics extracted from the time course of the membrane potential recording: the amplitude of the spike and the interspike interval. By using these quantities, the proposed method enables the accurate separation of excitatory and inhibitory contributions, yielding highly favorable results in the spiking regime. Author summaryQuantifying the continuous stream of inputs a neuron receives is key to understanding brain connectivity. Inside the brain, individual cells must maintain a tight balance between excitation and inhibition (E/I) to process information correctly, as any disruption in this equilibrium can impair its functionality. However, directly measuring the underlying excitatory and inhibitory synaptic conductances is technically challenging, and existing mathematical tools often fail when neurons enter their active firing regime. In this work, we introduce a novel computational strategy designed to extract and separate these time-varying conductances directly from the neurons spiking activity. By dynamically tracking just two accessible metrics - the amplitude of the spikes and the time intervals between them - our algorithm estimates both conductance profiles with high precision. Furthermore, we demonstrate that this procedure is highly robust against realistic experimental noise and data variability, providing an accessible framework that does not require complex hardware or an unfeasible number of repetitive experimental trials. By tracking changes in the E/I ratio of the synaptic input, this method provides an efficient approach to detecting pathological imbalances and understanding how local connectivity shapes cellular functionality.

Structural MRI is routinely acquired in clinical practice, yet quantitative morphometry has had limited impact on clinical decision-making. Overlapping symptoms, trajectories and comorbidities remain difficult to interpret within disease-specific frameworks, leaving it unclear how individual patients relate to the broader organization of brain disease. Here, we construct a cross-disease morphological reference space from 110,591 T1w MRI scans of 78,794 participants, spanning four disease families, 19 diagnoses, and seven subtypes. To construct this space, we developed NeuroMorph, an AI framework deriving thirteen interpretable morphological descriptors and individual normative deviation profiles. The reference space reveals shared and distinct morphological signatures that distinguish conditions within a hierarchy of disease families, diagnoses, subtypes and individual profiles. It identifies overlapping and comorbid morphological profiles and captures longitudinal deviations that precede clinical diagnosis and track progression. Together, these findings establish a unified framework for mapping brain disease organization and positioning individual patients within its morphological landscape.

Background: Acute cortisol responses to psychosocial stress vary substantially across individuals, yet how this variability is reflected in post-stress resting-state functional connectivity (rsFC) remains unclear. Although prior work has linked stress-related endocrine responses to brain connectivity, studies have been limited by small samples, region-of-interest approaches, or a sole focus on group-level analyses. Here, we investigated whether acute cortisol increase is associated with, and can be predicted from, whole-brain post-stress rsFC. Methods: We analyzed 339 healthy participants from two ScanSTRESS datasets using complementary inferential and predictive approaches. First, we used the Network-Based Statistic (NBS) to identify connected rsFC networks associated with acute cortisol increase, controlling for age, site, and sex/hormonal status. Second, we predicted participants' acute cortisol increase from their connectivity patterns using NBS-Predict and Connectome-Based Predictive Modeling (CPM). Together, we examined the cortisol-rsFC relationship at the population and individual levels. Results: Greater cortisol responses were associated with lower post-stress rsFC within a significant distributed network comprising 258 connections among 78 regions, centered on thalamic nuclei and pallidal regions and extending to default-mode, limbic, orbitofrontal, and cerebellar regions. Sex-stratified analyses revealed a significant negative association only in females, but formal sex-difference contrasts were not significant. NBS-Predict and CPM yielded modest but significant out-of-sample prediction, with predictive networks converging on subcortical and posterior cingulate regions. Conclusions: Post-stress rsFC carries convergent inferential and predictive information about individual HPA-axis reactivity. Stronger cortisol responses were characterized by reduced connectivity within a distributed subcortical-cingulate network, supporting a network-level perspective on neural-endocrine coupling following acute stress.

Action-selection is determined by a combination of goal-directed and habitual processes. Habits are defined as the reward-independent, stimulus-response relationships which form when an action is regularly executed in the same context, regardless of outcome. An influential computational model proposes that habit formation is driven by action prediction errors which occur when non-habitual actions are taken. It has been further suggested that action prediction errors are encoded in activity of specific dopamine neurons, and it has been recently observed that dopamine activity in the tail of the striatum follows a pattern consistent with the action prediction errors. However, the original models capture changes in habits across trials, but do not describe the time-course of action prediction errors within trials, hence it is difficult to directly compare them with dopamine activity. We begin by outlining the temporal-difference action learning algorithm, which uses biologically-plausible mechanisms to determine how dynamic changes in action intensity influence the resultant prediction errors across near-continuous time. We then demonstrate that dopaminergic data recently collected from the tail of the striatum is better represented by action prediction errors than reward prediction errors. Overall, our results support the existence of value-free action prediction errors and associated habitual behaviour in dopaminergic signals. Author summaryWhenever we choose one action over another, there are two ways that the selection can be made. We could take the time to consider what we want to achieve, calculate which action is the most likely to give us that outcome and balance it against the possible negative consequences. These goal-directed calculations are very time-consuming and our brains could not possibly do it for every choice. Instead, we often rely on the second method, habits, which learn to copy the actions that were most often chosen in the past. In this paper, we present a new model of learning that is based on biologically plausible brain networks and applies action prediction errors to update our habits across continuous time. Using simulations, we reveal testable predictions that are specific to our temporal-difference action learning model and build an intuition for its behaviour. Finally, this model is tested against real dopaminergic data from the tail of the striatum, and we show that our model provides better explanation for these data, than classic reward-based reinforcement learning models.

Oral contraceptives are used by millions of women worldwide, yet their cumulative effects on the female brain remain poorly understood. We analyzed resting-state fMRI from 192 women (never, current, and past users) using brain-dynamics metrics that quantify information transfer across spatiotemporal scales. Duration of use was associated with a progressive, age-independent modulation of brain dynamics: prolonged use was related to restricted local synchronization but enhanced information flow across scales. Machine learning predicted individual duration of use from brain dynamics, and node-level classifiers distinguished users from never-users based on a spatially distributed cortical pattern spanning multiple functional networks, most clearly detectable in long-term users. In past users, these associations were inverted in sign, scaling with prior duration rather than returning to the never-user baseline, suggesting active reorganization following cessation. Together, these results indicate that cumulative oral contraceptive exposure acts as a neuromodulator of information processing.